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As  the  complete  equations  for  plastic  flow  with  arbitrarily  large 
deformations  are  not  readily  available,  the  results  of  a  number  of  scattered, 
relevant  analyses  have  been  assembled  in  the  following  pages  to  provide  a 

summarj  of  the  appropriate  equations.  For  simplification,  tensor  notation 
is  used  everywhere. 


THE  CONSERVATION  LAWS 

The  equations  of  motion  for  a  fluid  are  given  below  for  a  general 
tensor  law  of  resistance: 

+  p  u^^^  =  0,  ntinuity,  (l) 


P 


momentum,  (2) 


P  ^  I  energy,  (3) 

where  the  convective  derivative  is  given  by  the  expression 

i  s 

Dt  "  dt  ^  ^i 

and  repeated  indices  are  summed.  Differentiation  with  respect  to  a 
coordinate,  x^,  is  denoted  by  preceding  the  index  with  a  comma,  e.g.. 


The  energy  equation  is  simplified  by  elimination  of  the  kinetic 
energy  term  in  Eq.  (3)  by  means  of  the  momentum  equation,  resulting  in 


P 


Dt 


ih) 


The  velocity  gradient,  u  ,  can  be  written  as  the  sum  of  a 

J 

symmetric  and  an  antisymmetric  part.  The  symmetric  nart,  e  ,  is  the 
strain-rate  tensor  and  the  antisymmetric  part,  is  the  angular 

velocity  vector  in  tensor  form. 


1 


2 


Thus^ 


where,  by  definition. 


u.  .  =  e.  .  +  u). 

i.]  ij' 


^  ^  ■‘“ji 

The  stress  tensor  is  always  s^Tnmetric,  a  result  found  by  computing  the 
total  moment  on  an  element  of  fluid  and  noting  that  the  surface  forces 
must  cancel.  Then, 

tj  “  “ji-  (T) 

Interchanging  dummy  indices,  using  Eqs.  (5)  and  (6),  and  noting  that  an 
expression  equal  to  its  negative  is  zero, 

“ij  "ij '  ti  “ji  * 

The  energy  equaulon  (i4-)  then  becomes 

»  i '  hj  hy  (8) 

The  stress  and  strain  tensors  can  be  written  as  the  sum  of  a 
hydrostatic  part  and  a  deviator  part,  defined  as  the  total  tensor  less 
the  hydrostatic  part.  The  total  tensors  have  been  denoted  above  by 
Greek  letters;  the  deviators  will  be  denoted  by  Roman  letters.  Then  the 
stress  tensor  becomes 


where 


‘’ij  “  -J’  hj  *  hj’ 


P  -  3  -  0, 


is  the  Kronecker  delta.  Similarly,  the  strain  rate  tensor  becomes 

'ij  '  3  «  hj  *  >  (11) 


where 


e  -  6^^,  .  0. 
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Using  thpse  relations^  the  five  flow  eq^uations  can  be  written  in 
terms  of  the  deviator  quantities: 


it  pe  = 


(13) 


Du. 

P  Dt“  =  ”P  S,.  .  ,  , 


DE 

=S,,  e,.. 


>1 

ij  '"ij 


W 

(15) 


THE  CONSTnUTIVE  EQUATIONS 

A  wide  variety  of  material  properties  can  be  described  by  a 
cons  itutive  equation  which  expresses  the  total  strain  rate  as  the  sue 
of  an  elastic  and  a  plastic  part: 


ij 

( total ) 


-i-  2 
2^1  ^ij 


b  3 


(.16) 


iJ  ^ 

(elastic)  (plastic) 

where  p  is  the  shear  modulus  of  elasticity;  b  is,  in  general,  a  function 

of  stress  and  strain.  Several  special  cases  in  which  the  second  stress 
invariant 

^2  '  (IT) 

appears  repeatedly  are  listed  below.  Where  k  appears,  it  denotes  the 
yield  strength  of  the  material. 

Classical  viscosit.y 

Shear  .modulus  =  p  =  <d;  ^  =  l/2  p,  p  =  dynamic  viscosity;  (18) 

Prandtl-Reuss  plasticity  equations 


)  (19) 


for  Jg  =  k^, 

0 

II 

.a 

for  Jg  <  k^; 

k 


Perznya  plasticity  equation^ 
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1  "  1 

for  Jg  k^, 

> 

b  =  0 

2 

for  <  k  ; 

/ 


Class  leal  elasticity 

b  =  0; 

Rlffid-plastlc  equations 

^  =  03,  vith  Eq.  (19)  or  Eq.  (20). 


(20) 


(21) 


(22) 


MATERIAL  ROTATION 

The  general  constitutive  equation,  (16),  is  only  valid  for 
coordinates  fixed  in  the  material.  If  it  is  to  be  used  in  connection 
with  a  set  of  Eulerian  flow  equations,  for  which  the  coordinate  axes  are 
fixed  in  space,  a  correction  for  rotation  of  the  material  must  be 
included.  The  physical  basis  for  this  correction  term  is  that  when  the 
material  rotates  there  is  a  change  in  the  con5)onents  of  the  stress  tensor, 
even  though  the  physical  state  of  stress  may  not  change.  Let  T. .  denote 
the  components  of  any  tensor  in  material  axes  and  T^^  the  components  in 
fixed  eoces .  Iben 


T.  .  =  T  a.,  a.  . 
ij  kX  ik  ji 

T  “  T  a,  .  a  .  . 
iJ  kX  ki  Xj 


(23) 

(24) 


and  the  direction  cosines,  which  are  time -dependent,  satisfy  the  identities 

(25) 

(26) 

Equation  (16)  ecomes,  in  this  notation, 

hj  ■  k  ®ij "  ^  ®ij  • 


ik  jk" 

a . .  a,  .  =  6 . 
jx  ki  jk 


tr 


It  is  required  to  replace  the  first  term  on  the  right  with  appropriate 

fixed-axis  derivatives.  Differentiating  Eq.  (2h) 

* 

^ki  ^  ^ki  '^kA  ^ki  '’XJ 

and  multiplying  by  a  .  a  . , 

pi  qj 

■"ij  ^pi  ~  "^pq  ^  \q  ^ki  ^pi  ^qj  j  ’ 


(28) 


The  quantities  -a^^^  a^^^  are  the  components  of  the  angular  velocity 
tensor,  .  To  see  this,  we  make  use  of  the  angular  velocity  vector, 
Q,  about  which  the  material  axes  are  rotating.  Let  a  denote  the  unit 
vector  along  the  A  axis  in  the  rotating  material.  The  angular  velocity 
of  the  endpoint  of  this  vector  is  in  vector  notation 

a  =  QX 


or  in  tensor  notation 


‘"ri  ®rjk 


(29) 


where 


IJiC 


I'Jj 

>  according  as  i,  j,  k  ^ 

j  form  an  even  I 
form  an  odd  ) 

1  o| 

do  not  form  a  { 

permutation 
of  1,2,3. 


(30) 


Multiplying  Eq.  (29)  by  a^^  and  using  Eq.  (26), 


\je  ^rji 


The  angular  velocity  vector,  Q  ,  is  related  to  the  angular  velocity  tensor 


ty 


u).r.  =  e  .  n, . 

ik  ijk  J 


Combining  these  two  equations  and  changing  indices, 

a  .  a . ,  =  -ID  , 
rl  i£  ri 

and  Eq.  (28)  beconKs 

T.  ,a,a.=T  +T,  +T  m  . 

1.3  pi  qj  pq  kq  ^  p/  '"/q 


(31) 


(32) 


(33) 
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^tultiplying  Eq.  (27)  by  eliminating  the  barred  quantities  by 

means  of  Eq.  (33)j  and  putting  the  deviator  stress  tensor,  S  ,  for  the 
general  tensor,  T.  ,  we  get,  finally, 

1 J 


ij 


®kj  \i "  ®ik  \i'>  * " 


{3^) 


The  term  on  the  right  of  the  energy-'  equation,  (15),  can  be 
obtained  by  multiplying  Eq.  (3^)  l>y  and  using  the  expression  for 
Jg  in  Eq.  (17).  The  terms  involving  m.  .  disappear  by  syns  Jtry 
considerations,  as  outlined  in  the  equations  b-elow: 


S.  u)  S.  , 
ir  rj  ij 


-  S .  f.i) .  S  .  . 
ir  jr  1.1 


-  S, .  u)  .  S.  =  0, 
ij  rj  ir 


The  Lagrangian  time  derivative  in  Eq.  (3^)  should  be  replaced  by  the 
convective  derivative  in  the  Eulerian  formulation;  the  energy,'-  equation, 
in  its  final  form,  then  becomes 


P 


Dt 


DJ. 


M  Dt 


+  b  J, 


2' 


(3^) 


The  first  term  on  the  right  is  the  rate  of  increase  of  distortional  strain 
energy  and  the  second  tern  is  the  work  done  by  the  plastic  strains.  The 
strain  energy  is  recoverable,  but  the  work  done  by  plastic  strains  is  not. 


THE  COIgLETE  EQUATIOI^ 

For  convenience,  the  equations  for  plastic  flow  are  assembled  below 
and  include  the  equation  of  state,  which  has  not  been  introduced  before. 

Continuity  Equation 


D£ 

Dt 


+  p  9  =  Oj 


(36) 


Momentum  Equation 


Du. 

"  dT  =  ■  "  Sij.j- 


(37) 


Energy  Equation 


(38) 
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^uation  of  State 


Constitutive  Equations 


P  -  F  (E,p) 


where 


®ij  2^  |~Dt’^  ^  '\±  ^Ik  \i 

®ii  - 


b  =  0 


Defining  Equations 


'  5  "  “j,i> 

"  2  ■  “j,l* 


i  e  6. . 

3  iJ 


for  =  k  , 


for  Jg  <  k 


(deviator  strain  rate), 

(^3) 

(rotation  rate),  (44) 


*  .  .  • 


THE  PLASTIC -RIGID  EQUATIONS 


This  situation  is  considerably  simplified  if  the  plastic-rigid 
assumption  is  introduced,  which  is,  in  effect,  to  neglect  the  elastic 
regime.  Then  the  constitutive  equation,  (4o),  becomes 


and  Eq.  (42)  is  sinqjly 


. . .  . .  -  ■»:'  . 11 
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Diff<=“rentiating  Eq,.  (^7)>  put+lng  ^  ~  P  (denoting  a  viscositylike  terra)^ 
and  making  some  straightforward  manipulation.. 


®ij,j 


P  (V^  u. 


+ 


.)  +  2  e 

i' 


ij 


m 


The  equations  are  similar  to  the  Navier-Stokes  equations,  hut  the 
"viscosity,"  n,  is  not  a  constant,  as  in  the  usual  viscous  flow. 
Summarizing  these  results,  we  have  for  plastic-rigid  flow  six  equations 
in  the  six  unknowns,  p,  u^,  p,  and  E.  The  last  three  equations  for  the 
variables  p,  0  are  not  essential  and  are  not  cotinted  in  the 

stn. 


§1  +  P  9  =  0, 

(50) 

^i  2  1 

(51) 

P  ^  ®  ~  ^  ^  \/®ij  ^ij  ^ 

(52) 

P  =  F  (E,p), 

(53) 

_  k 

^  ‘  2  /i— i-  ^ 

W 

®ij  "  2  ■  I  ®  ^iy 

(55) 

®  =  ^i,i* 

(56) 

CYLINDRICAL  COORDINATES 

The  equations  for  axially  symmetric  flows  are  given  below,  u  and 
w  denote  the  radial  and  axial  velocity  components,  r  and  z  denote 
coordinates  in  the  radial  and  axial  directions.  Dots  denote  ordinary 
time  derivatives  and  subscripts  denote  differentiation  with  respect  to 
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the  iiidJcated  coordinate, 

p  +  up^  +  wp^  +  pe  =  0, 

p  +  uu^  +  wu^)  -  -p^  +  iT  +  i  I 

+  2^.1  Pr  +  20i3  P,- 

p(w-^UW  +  ww)=-p  +u(u  +-W  +  W  +ifl) 

^  r  z'  ^z  ^  rr  r  r  zz  3  ®z' 


+  2e,_  p  +  e  _  p.  , 

13  r  33 


p  (E  +  uE^  +  wE^)  +  pe  =  k  , 


^11  =  -  3  ^ 

'31  ■  '13  '  5  (“i  '  -  3  ' 


P  ss  1/  •  ^ 

33  z  3  " 


9=u  +  —  +  w  , 

r  r  z  ’ 


2  .  1 


2.2.  ,uv^ 


e  -  u!  +  i  (u  +  w  )*=  +  +  (^) 

J  IJ  r  2  '  z  r  z  'r' 


3  ^ 


SPHERICALLY  SYMMETRIC  FLOWS 


In  the  case  of  spherically  s.’nranetric  flows,  the  equations  are  the 
following,  in  which  u  denotes  radial  velocity  and  r  denotes  distance 
from  the  origin. 


P  ■*■  ^Pj.  +  p9  =  0  ^ 


2  u 


p  (u  +  uu^)  =  -p„  +  p  (u_  +  )  +  2e^^  {1^  , 


rr  r 


. . . . I 
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